Restoration of the magnetic hc/e -periodicity in unconventional superconductors 
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We consider the energy of the filled quasiparticle's Fermi sea of a macroscopic superconducting 
ring threaded by an /ie/2e-vortex, when the material of the ring is of an unconventional pairing 
symmetry. The energy relative to the one for the hc/e- vortex configuration is finite, positive, 
and inversely proportional to ring's inner radius. We argue that the existence of this energy in 
unconventional superconductors removes the commonly assumed degeneracy between the odd and 
the even vortices, with the loss of the concomitant hc/2e periodicity in external magnetic field as 
a consequence. This macroscopic quantum effect should be observable in nanosized unconventional 
superconductors with a small phase stiffness, such as deeply underdoped YBCO with T c < 5K. 
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A fundamental property of all known superconductors 
is that their electrons form Cooper pairs. A direct mani- 
festation of this phenomenon is the quantization of mag- 
netic flux in units of hc/2e in multiply connected geome- 
tries [H . Such flux quantization may be considered as an 
effective spectroscopy of the charge of the carriers, and 
is often used as a proof of paired nature of the super- 
conducting state. A closely related phenomenon is the 
periodicity of various properties of multiply connected 
superconductors in the external magnetic field with the 
period that corresponds to the half flux quantum 0] . 

In this Letter we argue that this periodicity is in prin- 
ciple not exact in superconductors with unconventional 
pairing symmetry that support quasiparticle excitations 
at arbitrarily low energies. Our argument is qualitative 
and fundamental in nature, and based on a calculation 
of the difference in ground state energies of nodal quasi- 
particles of an unconventional superconductor in a pres- 
ence of a single (hc/2e) and a double (hc/e) vortex in 
an annular geometry. An estimate of this energy differ- 
ence indicates that the deviations from the familiar hc/2e 
periodicity may become directly observable in recently 
fabricated deeply underdoped cuprates. This manifesta- 
tion of macroscopic quantum coherence is a fundamental 
effect and raises the possibility of manipulation of spin 
currents (carried by quasiparticles) by controlled motion 
of magnetic fluxes (using, for example, a Hall bar) with 
applications to spintronics and other areas of applied sci- 
ence. 

The energy in question stems from the essential differ- 
ence in the interactions between the quasiparticles and 
hc/e- and ft.c/2e-vortices. Apart from the semiclassical 
Doppler shift of the quasiparticle energies common to 
both single and double vortices, the statistical, purely 
quantum Aharonov-Bohm (AB) phase of an hc/e- vortex 
can be exactly gauged away whereas the one of an hc/2e- 
vortex cannot. The ensuing topological frustration is 
felt by the quasiparticles arbitrary far from the cen- 
ter of the vortex via the AB gauge field which en- 
codes the sign change in the quasiparticle's wave-function 




FIG. 1: (Color online.) Annular system with the vortex in 
the nodal superconductor (gray region). 



as it is adiabatically dragged around the vortex. Vor- 
tices, their fluctuations and the concomitant AB and 
Doppler effects on quasiparticles in d-wave superconduc- 
tors have been a subject of much research in the past 
H @, 0, 0, S, 0, H El, E3- Here we consider the filled 
Fermi sea of nodal quasiparticles in an annular geometry 
(Fig. 1), and determine the excess in energy due to the 
AB gauge field of the ft.c/2e-vortex. We find a positive 
contribution to the condensation energy that derives pre- 
dominantly from the quasiparticles near the nodes and is 
inversely proportional to the hole radius R. For param- 
eters relevant to cuprates the excess energy is <~ 0.2K 
for a thick ring whose inner radius is about a microme- 
ter. Consequences for the quantization of magnetic flux 
in underdoped cuprates are briefly discussed. 

Let us assume a magnetic flux localized in the annu- 
lus made of an unconventional superconductor. At low- 
energies, the dynamics of quasiparticle excitations near a 
single node in the field of a vortex in the superconducting 
order parameter carrying a half flux quantum hc/2e may 
be described by the Hamiltonian 



H = v F (p x + a x )ai + v A {p y + a y )a 2 ~ mcr 3 , (1) 

where vf and va are characteristic velocities of the quasi- 
particle excitations in the two directions around a nodal 
point, en are the Pauli matrices. a(r) = (xe y —ye x )/(2r 2 ) 
is the AB vector potential resulting from the Franz- 
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Tesanovic Q (FT) gauge transformation in presence of 
the ft.c/2e-vortex in the order parameter. In |T]) we have 
included a gap m at the nodes for generality and set 
h = c = 1. We have also neglected the Doppler shift 
of quasiparticles, based on the following argument: the 
Doppler effect enters Hamiltonian ([I]) as another gauge 
field v(r) generated by the FT transformation. It always 
appears in the combination v(r) — (e/c)A(r), where A(r) 
is the electromagnetic vector potential. When the hole 
radius R (Fig. 1) grows to become comparable to the 
magnetic field penetration depth A, A(r) will screen out 
v(r), leaving the topological frustration encoded by a(r) 
as the sole long range effect. Thus, in practical situations, 
we expect the Doppler shift to be a secondary effect in 
rings of macroscopic size. 

We consider first the simpler case of isotropic velocities, 
vp = va- Setting the velocity to unity the eigenstates of 
the above Hamiltonian are found to be 



for I > 0, and 
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V g ( £-m)J,_ 1/2 (fcr)e^- 1 >'* 
iq y /q{E + m)J l+1/2 (kr)e il 1> 



(2) 
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y/g( E - m)J_i + 1/2 (fcr)e' 
-iq y /q(E + m)J_ l _ 1/2 {kr)e a 1> 



(3) 



for I < 0. Here the quantum number q = ±1 distin- 
guishes particle- like and hole-like states, I £ Z is angular 
momentum, k > is radial wave vector, and energy is 
E Qj k = qVk 2 + m 2 . Ji(x) are Bessel functions of the first 
kind. To find the above eigenstates it is necessary to 
regularize the gauge potential. Namely, the requirement 
of the square-integrability unambiguously determines all 
the eigenstates of the Hamiltonian |T]) , except ones with 
zero angular momentum. There are two I = states 
diverging as \j\fr at the origin, but keeping both of 
them leads to an overcomplete basis in the Hilbert space. 
On the other hand, requirement of non-divergence of the 
states at the origin is too restrictive, since it leads to an 
incomplete eigenbasis. Only linear combinations of the 
two states specified by a single parameter are allowed [l3| 
but in order to select a single eigenstate from all allowed 
states, the gauge potential has to be regularized. Here, 
we considered the vortex as a cylinder of a finite radius R 
in which the AB flux is uniformly distributed. By match- 
ing the solutions inside and outside the cylinder, and tak- 
ing the limit R — > 0, we found that the eigenstate in the 
zero angular momentum channel has the form given by 
©, with the lower component diverging at the origin. 
This is in agreement with the result of an alternative 
regularization 11, T3 |. 



Let us now calculate the local density of states (LDOS) 
for gapless nodal quasiparticles, m = 0, defined as 



p(e,r) J dk\y q , l<k \ 2 6(e-E q! , 

q-l 



(4) 



Using the eigenstates given by Eqs. |(2]) and (J3j) , we re- 
produce the LDOS of Ref. [13 in the form 



P(e, r ) 



cos(2|e|r) , |e| 



27r 2 r 



(5) 



The expression for the LDOS can further be simplified 
using the little known Mitrinovic identity 14 1 



£[J p+i (x)] 2 =p / -[J p (t)} 2 ~ Mx)} 2 - (6) 
1=1 • )q 

For p = 1/2, this yields the form of the LDOS which is 
more convenient for the later calculations 
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P(e,r) = -o 



cos(2|e|r) 
2r 



|e|Si(2|c|r) 



(7) 



where the standard sine-integral function is defined as 
Si(x) = L dt sint/t. In the vicinity of the vortex, in 
the region r <§; l/|e|, the LDOS diverges as 1/r. This 
behavior of the LDOS originates from the states in the 
zero angular momentum channel that diverge as l/^r, 
when r — > 0. On the other hand, far from the vortex, 
p(e,r) — > po(ej r ) == N/27]", as in the vortex-free system. 
Of course, the LDOS in the system with the vortex carry- 
ing an integer number of the flux quanta, nhc/e, n £ Z, 
is the same as in the free system. Namely, the vector 
potential corresponding to n/ic/e-vortex is 2na(r), and 
in that case the eigenstates of ([I]) have the form 
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y/q(E - m)J|,_ 1 +n | (kr)e^ 1 ^ 
tq^/q(E + m)J ll+nl (kr)e u * 



(8) 



The LDOS for gapless nodal quasiparticles po(e, r) is then 
uniform and independent of the integer n, as required by 
gauge invariance. 

Starting with the compact form of the LDOS, we may 
compute the energy cost of having an /ic/2e-vortex by 
integrating the Eq. J7]) over the energy and the area of 
the ring in Fig. 1. This procedure should be accurate 
for a ring of a macroscopic size, when the effects of the 
boundaries and of the discreetness of the spectrum be- 
come negligible. The DOS for the ring is then 

p(e) = J d 2 rp(e, r) = I(R + I, e) - I(R, e), (9) 
cos(2ei?) sin(2ei?) 



where 
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8i(2eR) 



2eR (2eRf 



(10) 
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R + I and R are the radii of the outer and the inner 
annulus, respectively. The total energy of the system 
then becomes 



£ 



deep(e)=£(R + l)-£{R), (11) 



with 



S(R) = - 



A 3 R 2 



Si(2i?A) 



cos(2i?A) sin(2i?A) 



3tt 

1 - cos(2i?A) 
4i? 3 A 3 



2i?A 



(2RAy 



(12) 



Here, A is a high-energy cutoff, and the minus sign in 
Eq. (fTT]l takes into account that only hole-like states are 
occupied in the ground state. Using the asymptotic form 
of the sine-integral function for large values of its argu- 
ment, we find the energy cost of an ft.c/2e-vortex in the 
annulus of a macroscopic size R^> 1/ A to be 

f ^ f - fo =i2^To( 1 + ° ( ^)' (13) 

where the total energy for the /ic/e-vortex (or the vortex- 
free system) is £ = A 3 [i? 2 — (R + Z) 2 ]/6. When the 
thickness of the annulus is much larger than its inner 
radius, I » R, the extra energy cost due to the presence 
of an /ic/2e-vortex in the order parameter, to the leading 
order in 1/AR and R/l, is simply 



f 



hvp 
YhxR' 



(14) 



where we have also restored Planck's constant and the 
Fermi velocity vp — va previously set to one. The energy 
cost for having an /ic/2e-vortex in the system therefore is 
positive, and in the macroscopic limit, for a thick annulus, 
inversely proportional to its inner radius. Notice that 
the leading term in £ v is independent of the high-energy 
cutoff A, in accord with our assumption that the effect is 
due to the low-energy quasiparticles near the nodes. The 
long-wavelength, linearized, description we postulated in 
Eq. (1) is thus internally consistent for an annulus of a 
macroscopic size. The result in Eqs. lfl3"l) and (fl4|l also 
reflects the fact that the presence of the vortex affects 
the LDOS in its (macroscopic) vicinity the most. 

We can now turn to the general and a physically more 
relevant case when the two characteristic velocities of the 
nodal quasiparticles, vp and va, are different. By rescal- 
ing the coordinates, x' — x/vp, y' — u/va, and choosing 
a gauge such that the vector potential has the same form 
as below Eq. (1) in the new coordinates (x',y'), Hamil- 
tonian |l|) may be transformed to a form with isotropic 
velocities 0, [l(3, EH. The rescaled momenta are now 
k' x = vpk x , k' y = VAky, and the dispersion assumes an 

isotropic form, E q ^ = qk', with k' = J k' x 2 + k' y 2 . The 
LDOS then becomes p(e,r')/2w, with the extra factor of 




Va A 




V F A 



k: 



FIG. 2: Shape of the Brillouin zone before and after rescaling 
of the momenta. 



2ir arising from the elliptic shape of the Brillouin zone in 
the new coordinates, shown in Fig. [2J 

The total energy of the anisotropic vortex-free system 
in the annular geometry is then £q — G(R + I) — G(R), 
where 



G(R) = - [ d 2 r [ d 2 k'p (k' 

Jr<R 7s(A) 



(15) 



E(A) is the Brillouin zone in Fig. 2, and po(k') = 
k' /(2ir) 2 is the LDOS of the vortex-free system. The 
integration yields the total energy of the flux-free system 



,31 vi 



1 3 
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VF > VA, 



1;1 



7 A 



(16) 



and 2F\{a,b;c;x) is the hypergeometric function. If 
va > vp, the two velocities should be exchanged. In 
the presence of an /ic/2e-vortex the total energy of the 
system is thus 



£ = g{R + i)-g{R), 



where 



G[R) 



d 2 r' 



Cl(R) 



d 2 k'k'p(k',r') 



2(A) 



(17) 



(18) 



v A y' 2 < R 2 - To the leading 



Q(R) is the ellipse v 2 F x' 2 
order in 1/ ( AR) , one then finds £ = £q + £ v , where £q is 
the energy of the flux-free system given by Eq. fTS)! , and 
the contribution to the total energy arising solely from 
the presence of the AB vector potential when vf > va is 



£ v = 



hlvF 



l2irR(R + l) 



11 v 

r-2^-v 



(19) 



In the isotropic case, vf — va — 1, we obtain the re- 
sult Cd, since a-PiQ-A -1/2; 1;0) = 1. In the opposite 
limit of a large velocity anisotropy, the energy cost is de- 
termined by a larger of the two velocities, because the 
function 2-Fi(l/2, — 1/2; l\x) is monotonic and bounded 
on the interval [0,1], 2/tt = 2^(1/2, -1/2; 1; 1) < 
2*1(1/2, -1/2; l;x) < 2*1(1/2, -1/2; 1; 0) = 1. 
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The result in Eq. (fl"9| pertains to the energy of filled 
quasiparticle Fermi sea. However, we can straightfor- 
wardly import it into the fully self-consistent computa- 
tion of the total superconducting condensation energy: 

|A| 2 

E tot (A, hc/2e) = E qp (A, hc/2e) + J— L , (20) 

where E qp (A, hc/2e) is the energy of the Fermi see with 
an hc/2e vortex and gap parameter A and g is the ef- 
fective coupling constant. We have assumed that A is 
essentially uniform since any non-uniformity enters only 
on the microscopic scale <C R or I. By adding and sub- 
tracting E qp (A, 0) we obtain: 

E tot {A,hc/2e) = £ v + E tot {A,0) , (21) 

where £ v is given by Eq. (fl9| . By minimizing 
E to t{A, hc/2e) with respect to A we obtain the total con- 
densation energy in presence of hc/2e vortex. This gives 
A = A + 8 A, where A is the value that minimizes 
iJ tot (A,0), and 5 A oc £ v . For macroscopic R, I » 1/A, 
£ v is arbitrarily smaller than E to t{A, 0) and SA <C A . 
This implies that, to the leading order, the presence of an 
hc/2e vortex increases the condensation energy by pre- 
cisely £ v l|19p. with A = Ao; the leading correction is 

For a finite s-wave gap vnR <g; 1 the calculation is simi- 
lar but considerably more cumbersome. We find that the 
excess energy in Eq. lfT4|) decreases with the gap m, and 
is essentially zero already for mR ss 1, which, crudely, 
would correspond to an inner radius of a micrometer in an 
aluminum ring. This is in accord with our interpretation 
of the effect as being due to the nodal quasiparticles, and 
with the high accuracy of the observed /ic/2e-periodicity 
in standard low-T c superconductors. 

We can estimate the above energy cost of the 
/ic/2e-vortex from the values of the Fermi velocity 
and the velocity anisotropy in YBa2Cu307_a (YBCO) 
and BisSrsCaCuaOs+y (BSCCO) obtained by ARPES 
[lB | and the thermal conductivity measurements 
ARPES yields a value of the Fermi velocity vp ~ 3 • 
10 5 m/s, which appears to be universal in cuprates. The 
velocity anisotropy is vf/va ~ 14 in YBCO, while in 
BSCCO vf/va ~ 19, yielding the total energy penalty 
of having an /ic/2e-vortex, £' ot = N£ v ~ 0.2K, for an 
annulus with the inner radius R = 1/im, and N = 4 as 
the number of nodes in a d-wave superconductor. 

The finite energy cost of an hc/2e- vortex will affect 
the quantization of the magnetic flux when it becomes 
comparable to the second relevant energy scale in the 
problem, namely the superfluid density, p(T). By lift- 
ing the parabolas centered at the hc/2e flux in the text- 
book energy vs. magnetic flux plot [l[ by £ v it is easy 
to see that the width of the hc/e relative to the one of 
the hc/2e plateau becomes longer by an amount propor- 
tional to 8 = £ v /p(T). This is typically a small number: 



in optimally doped YBCO, for example, S w 10~ 4 . Re- 
cently, however, single crystals [ijj and thin films of 
severely underdoped YBCO have been studied with the 
unprecedented low p(0) ~ IK, when expressed in energy 
units This extremely underdoped regime where the 
phase stiffness can be rendered arbitrary small with un- 
derdoping offers the best chance for an observation of the 
asymmetry between even- and odd-flux vortices predicted 
in this paper. 

The asymmetry between hc/e- and hc/2e- vortices in 
unconventional superconductor has also been recently 
found in the numerical solution of the Bogoliubov - de 
Gennes equations for a mesoscopic superconducting loop 
[13] (see also (Hi)- The results and conclusions of this 
work seem broadly consistent with ours. 
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